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ABSTRACT 


The  rate  o£  evaporation  of  a  pure  liquid  at  tempe¬ 
rature  T  into  a  vacuum  is  shown  to  be  equal  to  a  *  c  .  n 

s 

where  a  is  a  coefficient  depending  on  equilibrium  molecular 
properties,  c  is  one-fourth  the  gas-kinetic  mean  velocity  of 
the  molecules  at  temperature  T  and  ng  is  the  number  density 
of  the  vapor  molecules  in  equilibrium  with  the  liquid  at  T. 

This  result  is  not  new.  Its  derivation  and  tinder- 
lying  assumptions  are  of  considerable  theoretical  interest. 
Previous  derivations  which  have  led  to  similar  or  slightly 


different  results  are  discussed 
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Introduction 

In  1832,  from  gas  kinetic  theory.  Hertz  (1)  arrived 
at  the  conclusion  that  there  must  exist  a  maximum  rate  of  evapo¬ 
ration  rg  from  a  liquid  surface  at  a  temperature  T  into  a 
vacuum.  According  to  Hertz: 

«  d  •  ng  (1) 

where  c‘  is  one-fourth  the  gas-kinetic  mean  molecular  speed  at 
temperature  T  and  ng  is  the  number  density  of  molecules  of 
vapor  in  equilibrium  with  the  liquid  at  T. 

This  result  of  Hertz  is  of  special  interest  in  that  it 
showed  for  the  first  time  how  a  rate  process  could  be  calculated 
from  gas  kinetic  theory  supplemented  by  a  fundamental  assumption 
of  equilibrium. 

Indeed,  the  reasoning  behind  equation  (1)  is  as  follows. 
Consider  a  liquid  in  equilibrium  with  its  vapor.  From  gas-kinetic 
theory,  the  frequency  of  collisions  of  vapor  molecules  with  the 
surface  is  given  by  c  .  n  .  If  every  collision  leads  to  con- 

O 

densation,  c  «  ng  represents  also  the  rate  of  condensation  and 
it  is  the  maximum  rate  of  condensation.  Since  equilibrium  pre¬ 
vails,  0  .  ng  is  also  the  maximum  rate  of  evaporation.  If  now 
all  molecules  in  the  vapor  phase  are  removed  so  that  the  liquid 
is  allowed  to  evaporate  freely  into  a  vacuum,  it  is  assumed  that 
the  maximum  rate  of  evaporation  is  still  given  by  the  equilibrium 

s* 


value  c  .  n 


-2- 


If,  on  the  other  hand,  as  first  noted  by  ICnudsen  (2), 
only  a  fraction  a  of  the  molecules  striking  the  surface  at 
equilibrium  actually  condenses ,  the  maximum  rate  of  condensation 
and  evaporation  at  equilibrium  is  a  •  c  «n  .  Again,  if  it  is 

5 

assumed  that  the  same  rate  obtains  away  from  equilibrium,  the 
maximum  rate  of  evaporation  into  a  vacuum  becomes: 

rQ  »  a  •  c  .  ng  (2) 

where  a  is  the  so-called  evaporation  or  condensation  coefficient. 
It  is  clear  that  a  must  have  the  same  value  for  condensation 
and  evaporation  since  it  was  first  introduced  to  describe  an 
equilibrium  process.  Then  by  virtue  of  microscopic  reversibility, 
a  must  be  the  same  in  both  directions.  On  the  other  hand,  if 
evaporation  ceased  to  be  an  equilibrium  process  away  from  equili¬ 
brium,  the  value  of  a  might  very  well  change  and  be  different 
for  condensation  or  evaporation. 

The  expression  of  Hertz  as  modified  by  Knudsen  was 
first  verified  by  the  latter  in  a  study  of  the  evaporation  of 
mercury.  Knudsen  (2)  found  a  ■  1  for  mercury.  This  result  was 
confirmed  by  Volmer  and  Estermann  (3)  who  found  a  *  1,00  +  0.07 
for  mercury  in  an  extended  temperature  range  where  the  rate  of 
evaporation  varied  by  a  factor  of  more  than  10  . 

Thus,  the  theoretical  result  of  Hertz  rests  on  solid 
experimental  grounds.  There  still  exists  no  other  case  where  a 


rate  process  has  been  calculated  and  the  calculated  value  sub¬ 
mitted  to  such  a  close  experimental  verification.  In  particular, 
the  fundamental  assumption  of  equilibrium  must  be  considered  as 
proved  by  the  careful  data  of  Volmer  and  Estermann. 

In  the  case  of  molecules  with  internal  degrees  of  free¬ 
dom,  gas-kinetic  theory  may  be  inadequate  and  it  is  of  interest 
to  apply  absolute  rate  theory  to  the  problem  of  evaporation.  The 
importance  of  rotational  motion  in  this  problem  was  first  pointed 
out  by  Herzfeld  (4).  In  what  follows,  the  rate  of  evaporation 
will  be  derived  from  the  usual  form  of  absolute  rate  theory.  This 
derivation  will  then  be  compared  to  other  derivations  in  the 
literature. 


Rate  of  Evaporation  at  Equilibrium 

The  simplest  approach  is  to  repeat  the  reasoning  of 
Hertz.  If,  at  first,  the  molecule  considered  has  no  internal 
degrees  of  freedom,  the  result  of  Hertz  must  of  course  be  repro¬ 
duced  . 


In  a  system  containing  liquid  and  vapor  in  equilibrium. 


the  rate  of  condensation  r  will  be  given  by  the  usual  formula 

C 


of  absolute  rate  theory: 
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RT 


(3) 
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Notations  are  as  follows: 

v  :  Eyring  *  a  frequency  :  lcT/h 

± 

F  :  the  partition  function  of  the  transition  state  per  unit 
interfacial  area 

F  :  the  partition  function  of  the  gas  molecules  per  unit 
g 

volume 

± 

Eq  :  the  activation  energy  at  absolute  zero. 

The  transmission  coefficient  has  been  assumed  to  be 
equal  to  unity.  If  there  is  no  activation  barrier,  (3)  reduces 
to: 


r 


c 


■  v 


(4) 


The  assumption  will  now  be  made  that  the  transition 

state  still  possesses  translational  mobility  in  two  dimensions, 

parallel  to  the  surface.  Since  the  molecules  have  no  internal 

degrees  of  freedom,  both  F  and  F  have  translational  contri- 

8 

butions  only,  respectively  in  two  and  three  dimensions.  Therefore: 


F*  -  (v/c)2  (5) 

F  *  (v/c)2  (6) 

O 

Substitution  of  (5)  and  (6)  into  (4)  gives: 

r  ■  c  .  n 

c  s 

But  since  the  liquid  is  in  equilibrium  with  the  vapor, 

r  »  r  .  If,  as  in  the  classical  theory,  it  is  now  assumed  that 
c  e 

this  is  also  the  rate  of  evaporation  when  all  molecules  in  the 
vapor  are  removed: 
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r  *  c  «  n 
e  s 

an  expression  identical  to  (1.)  as  expected.  The  case  of  interest 
is  that  where  the  molecules  have  internal  motion  characterized 
in  the  gas  and  in  the  transition  state  by  the  internal  partition 


functions  f  and 

S 

f* 

respectively. 

Then: 

F* 

an 

(v/c)2  f* 

(7) 

Fg 

m 

(v/c)3  fg 

(8) 

Hence,  the  rate  of  evaporation  becomes: 

re  *  (f+  /*g>  •  c  •  t.s  (9) 


A  final  assumption  is  to  take  for  the  internal  partition 
function  of  the  transition  state  the  value  of  the  internal  parti¬ 
tion  function  of  the  liquid,  f^. 

The  result  is  now: 

re  "  <fL/fg>  c  *  ns  (10) 


It  is  tempting  to  identify 
coefficient  a  : 


a  * 


fT/f  with  the  evaporation 
«  8 

(ID 


Calculation  of  a  is  then  feasible  in  principle  from 


equilibrium  properties  of  the  liquid  and  its  vapor. 
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Rate  of  Evaporation  into  a  Vacuum 

The  same  result  should  be  obtained  without  recourse  to 
the  circuitous  route  imagined  by  Hertz.  But  then  it  is  necessary 
to  use  a  particular  model  of  the  liquid  although  the  details  of 
this  model  must  disappear  in  the  final  expression  for  the  rate  of 
evaporation. 

Following  essentially  the  reasoning  of  Fenner  (5) ,  the 
calculation  of  the  rate  will  be  done  on  the  basis  of  the  free 
volume  theory  of  ICinkaid  and  Eyring  (6) . 

First,  the  relation  between  free  volume,  saturation 
pressure  and  heat  of  vaporization  will  be  recalled.  Then  this 
result  will  be  applied  to  the  calculation  of  the  rate  of  evapo¬ 
ration. 

The  complete  partition  function  for  one  mole  (N 

molecules)  of  the  liquid  is: 

QL  ”  If  {<v/c>3  Vf  exP  (E/RT)  *  fL  }  (12) 

where  is  the  free  volume  of  the  liquid,  per  mole  and  E  is 

the  molar  internal  energy  of  vaporization  at  absolute  zero. 

On  the  other  hand,  the  complete  partition  function  for 
one  mole  of  ideal  vapor  occupying  a  volume  V  is: 

Vtfr{<v/0)3v-  fg}N  (13) 

The  Helmoltz  free  energy  F  per  mole  is  related  to  the 
complete  partition  function  Q  : 
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F  *  -k  T  in  Q 


(14) 


The  molar  Gibbs  free  energy  or  chemical  potential  4  of  liquid 
and  vapor  are  then  given  by  the  expressions: 


“l  ■ 

-k  T  In  QL  +  p  VL 

**-k  T  in 

(15) 

"•«  ■ 

-k  T  in  Qg  +  p  V 

-  -k  T  in  Q„  4*  RT 
g 

(16) 

since  the  term  p  for  a  condensed  phase  can  be  neglected. 

At  equilibrium  between  liquid  and  vapor  °r 


V-  *  f.  •  exp(E/RT) 

in  ,  -  -  -1 

g 

Therefore,  putting  as  before  (f. /f_)  ■  a  : 

O 


a  •  -1  .  e'1  •  e-E/RT 


(17) 


(13) 


Introducing  volumes  per  molecule  v  «  V/N  and  v^  * 
V^/N  and  noting  that  v  *  (l/ng) »  we  have  finally: 


v 


f 


1  -E/RT 

1,1  .  © 

a  •  e  •  n 


(19) 


s 

To  calculate  the  rate  of  evaporation,  the  following 
usual  assumptions  are  now  made:  the  transition  state  is  in  equi¬ 
librium  with  reactants  and  the  transmission  coefficient  is  unity. 
Then: 


* 

r  »  v  •  n 

e 


(20) 


where  n*  is  the  number  of  activated  complexes  per  unit  inter¬ 
facial  area.  The  activated  complex  is  assumed,  as  previously,  to 


be  free  to  translate  in  two  dimensions  parallel  to  the  surface. 
In  the  calculation  of  n^,  the  translation  along  the  reaction 
coordinate  must  be  omitted.  Finally,  the  internal  partition 
function  of  the  activated  complex  is  assumed  to  be  identical  to 
that  in  the  liquid, 

di 

Then  the  complete  partition  function  Q  for  the  tran¬ 
sition  state  is: 

Q t  “  -|(v/c)2  .  f L  .  a  j-  (21) 

where  a  is  the  area  available  to  each  activated  complex  and 
therefore  a  *  l/r& .  The  equilibrium  condition: 


where  u  is  the  chemical  potential  of  activated  complexes  gives 
now: 

(v/c)  V£  *  exp  (E/RT) 

in - -  "  0  (22) 


Hence : 


(v/c) 


(23) 


since  by  Stirling's  approximation: 
N*  -  (N/e)N 


(24) 
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Finally: 


±  c  -E/RT 

r  -  v  .  vl  -  - — -  e 

e  e  .  v- 


(26) 


The  free  volume  per  molecule  in  the  liquid  can  now  be 
eliminated  by  substitution  of  its  value  derived  above  (19) : 

c  -E/RT  °  *  e  •  n. 


re  mt  e 


-e/rt 


(27) 


As  expected,  the  result  is  identical  to  that  obtained 
previously  (10) : 


r  »  a  *c  *  n 
e  s 


(23) 


Discussion 

In  the  final  result,  equation  (28) ,  any  reference  to 
the  particular  model  of  the  liquid  has  disappeared  and  the  ex¬ 
pression  for  the  rate  of  evaporation  is  the  same  as  equation  (10) 
which  was  derived  following  the  method  of  Hertz  without  any  assump¬ 
tion  as  to  the  nature  of  the  liquid  state. 

If  a  different  assumption  is  made  concerning  the  modes 
of  motion  of  the  liquid,  the  final  result,  equation  (28)  will  not 
be  affected  because  its  derivation  rests  on  the  two  equilibrium 
relations : 


"L  ■  •‘g 


between  which  the  value  of  the  free  volume  is  eliminated. 
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For  example,  in  the  treatment  just  given,  it  has  been 
assumed  implicitly  that  the  liquid  possessed  its  full  communal 
entropy.  Indeed,  the  partition  function  for  the  liquid  given  by 
equation  (12)  gives  to  each  molecule  the  possibility  of  roaming 
through  the  total  free  volume  V^. 

If  it  is  now  assumed,  for  instance,  that  each  molecule 
is  restricted  to  motion  within  its  own  cell,  the  communal  entropy 
vanishes  and  (12)  becomes: 

%m{  (v/c)3  vf  exp  (E/RT)  .  fL  ^  (12') 

Then  the  condition  for  equilibrium  between  liquid  and 
vapor  gives  for  the  free  volume  per  molecule  a  value  e  times 
larger  than  previously: 


a  •  n. 


e 


-E/RT 


(19*) 


But,  correspondingly,  the  condition  for  equilibrium 


between  liquid  and  transition  state  will  lead  to  a  value  for  n 
which  is  also  e  times  larger  than  previously: 


+ 


1 

vf  (v/c) . 


-E/RT 

e 


(25*) 


so  that  the  final  result  (28)  is  the  same  as  before.  However, 
Penner  in  his  original  derivation  (5)  used  essentially  equation 
(25')  together  with  equation  (19):  this  led  him  to  a  value  for  r 

e 


rfi  -  e  .  a  .  c  .  ng 
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which  is  of  course  e  times  too  high.  He  then  speculated  that 
this  value  should  be  reduced  by  a  transmission  coefficient  smaller 
than  unity,  expressing  the  possible  failure  of  the  equilibrium 
hypothesis.  The  situation  has  been  clarified  (7)  and  it  is  now 
clear  that  the  correct  expression  (28)  can  be  obtained  without 
recourse  to  any  ad  hoc  assumption  concerning  the  fundamental 
equilibrium  hypothesis  or  the  transmission  coefficient. 

More  recently,  Mortensen  and  Eyring  have  derived  once 
more  equation  (28) :  their  method  is  somewhat  confusing  in  our 
opinion  because  they  identify  a  *  fL /fQ  with  the  transmission 
coefficient  in  the  customary  rate  equation.  This,  of  course,  is 
not  necessary  at  all.  In  particular,  as  shown  above,  evaporation 
can  be  treated  like  a  perfectly  normal  rate  process  with  a  trans¬ 
mission  coefficient  equal  to  unity  and  the  usual  equilibrium  be¬ 
tween  activated  state  and  reactants. 

Since  a  can  be  evaluated  in  several  ways,  as  shown 
by  Mortensen  and  Eyring  (8) ,  it  becomes  highly  interesting,  from 
the  standpoint  of  theoretical  chemical  kinetics,  to  accumulate 
reliable  experimental  data  on  the  rate  of  evaporation  of  liquids. 

This  has  been  undertaken  in  our  Laboratory  and  results 
obtained  with  glycerol  have  been  discussed  elsewhere  (9) .  Further 
work  is  now  in  progress,  using  other  associated  liquids  for  which 
values  of  a  smaller  than  unity  are  expected  theoretically.  The 
possibility,  in  such  studies,  of  throwing  more  light  on  the  nature 
of  the  transition  state,  appears  quite  attractive. 
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